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Abstract 

T^j-" ' Operator method and cumulant expansion are used for nonperturbative calculation of the partition 

| function and the free energy in quantum statistics. It is shown for Boltzmann diatomic molecular 

, gas with some model intermolecular potentials that the zeroth order approximation of the proposed 

£\j ■ method interpolates the thermodynamic values with rather good accuracy in the entire range of 

both the Hamiltonian parameters and temperature. The systematic procedure for calculation of the 
corrections to the zeroth order approximation is also considered. 
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I. INTRODUCTION 

At present, most physical problems of interest can be solved neither exactly nor in terms of the canonical pertur- 
bation theory (CPT), therefore a considerable attention is given to the development of methods of nonperturbative 
analysis of quantum systems (see, e.g., review 1] and references therein). One of such methods, called the operator 
method (OM) of approximate solution of Schrodinger equation, proved to be very effective. The OM waspresented in 
and approved for a number of physical problems including those with many degrees of freedom 3-0. The main 
advantages of the OM are defined by the fact that its zeroth order approximation allows one to calculate the Hamil- 
tonian eigenvalues and eigenvectors with rather high accuracy in the entire range of both the Hamiltonian parameters 
and quantum numbers of the states. Besides, the further approximations of the OM lead to the sequences converging 
to the exact solutions. However, most of the particular applications of the OM have been considered before for the 
systems in "pure" quantum states. 

The nonperturbative methods are of great interest also in quantum statistics, when calculating thermodynamic 
characteristics of quantum systems. These characteristics can be expressed through either the partition function or 
the free energy of a system. Most of known nonperturbative methods in quantum statistics are based Feynman path 
integrals (e.g., [tj and references therein). However, an alternative representation of the partition function based on 
the direct summation over the energy states is also important because it could be more convenient and descriptive 
for some applications, especially in the atomic and molecular physics when the system is characterized by a large but 
finite number of degrees of freedom. For example, the additional procedure of bosonization is necessary when the 
functional integrals are used for the spin systems, and that leads to complication of the real Hamiltonian (see, e.g., 

The purpose of this paper is to generalize the OM for the nonperturbative calculation of the thermodynamic values 
in quantum statistics with the Schrodinger representation for quantum systems. A specific feature of this problem 
in comparison with application of the OM for "pure" quantum states is that the partition function is defined by an 
energy spectrum as a whole and include the summation over all states. As a result, the accuracy of calculation of 
the thermodynamic values for some temperature interval can be rather low even in the case when the Hamiltonian 
eigenvalues are found with high precision (see example in Sec. 2). Therefore we should analyze whether the OM 
accuracy for the spectral problem Q ensures the calculation of the thermodynamic values in the entire temperature 
range. Besides, some additional procedure should be developed for approximate summation over the states. It is shown 
in the paper that the regular nonperturbative method for calculation of the thermodynamic values can be developed 
on the basis of the OM combined with the cumulant expansion (CE) In order to approve our approach and 

compare it with other approximate methods, the Boltzmann diatomic molecular gas is considered and the partition 
function and the free energy corresponding to the molecular internal degrees of freedom arc calculated. It is shown for 
some particular examples that the zeroth order approximation of the proposed method leads to the uniformly suitable 
estimation for the thermodynamic values in the entire range of both temperature and Hamiltonian parameters. The 
systematic procedure for calculation of the consequent approximations is also considered. 

The paper is organized as follows: in Sec. 2 we discuss some general definitions of uniformly suitable estimation 
that concerns any nonperturbative approach. The procedure for nonperturbative calculation of the Hamiltonian 
eigenvalues on the basis of the operator method is described in Sec. 3 as the first part of the proposed approach. The 
approximate summation over the states by means of the cumulant expansion is considered in Sec. 4 as the second part 
of our method. In order to approve this expansion, the analytical formula for interpolation of the rotational partition 
function of the diatomic molecular gas is obtained in the same section on the basis of the method proposed. In Sec. 5 
we use our approach for calculation of the partition function and the free energy for an equilibrium system of the 
anharmonic oscillators and compare the calculated values with the results of other methods. The obtained results 
are of interest not only as the approval of the considered method, but also for the description of thermodynamic 
characteristics of real molecular systems when anharmonic effects should be taken into account. 

II. FORMULATION OF THE PROBLEM 

It is well known that if the Hamiltonian of some quantum system includes a small parameter, the regular procedure 
of calculation of the observable characteristics of this system can always be constructed in the form of a power series 
in this parameter. For example, it can be the canonical form of the perturbation theory with expansion in a power 
of small coupling constant (g <C 1); an expression as powers in g^ 1 for the strong coupling limit; a series in terms 
of powers of Ti in the quasi-classical approximation; low- or high-temperature expansions for thermodynamic values. 
However, in most cases these expansions are the asymptotic ones, and hence cannot be used directly in a wide range 
of the Hamiltonian parameters. On the contrary, the main goal of the nonperturbative methods is to calculate 
physical characteristics of a system in the entire range of its parameters. Further in the paper, when we refer to the 
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nonperturbative calculation of some physical values, we assume that either there is no small parameter in the system, 
or its value lies out of the domain of applicability of the asymptotic expansions - . 

It is convenient to define the concept of the "uniformly suitable estimation" (USE) in this case. Let us consider 
some physical variable with eigenvalues F n (X) depending on the quantum number n and the physical parameter A 
(the totality of the parameters and quantum numbers are implied for systems with many degrees of freedom) . Let us 
introduce 

Definition 1: the function F n °\X) is the USE for F n (X) if the following inequality holds in the entire range of 
variation of the values n and A 



Fn'{X) — F n (X) 

F n (X) 

Here the parameter < 1 is assumed to be independent of n and A, and defines the accuracy of the USE. We also 
consider 

Definition 2: there are the sequence of functions F^ S \X) and the method for their calculation corresponding to the 
decreasing sequence of the parameters ^ ; s = 0, 1, 2, . . ., so that 



(0) 



(1) 



lim FW(A) = F n (A). (2) 

s — >oo 

It seems that Definition 1 is not constructive because the exact values F n (X) are unknown in general case. However, 
there are several possibilities to estimate the value In particular, one can compare the asymptotic series for 

F n (A) in various limit cases either of the parameter A or of the quantum number n with the corresponding expansions 

of the function Fn (A) taking into account that inequality should hold for all cases simultaneously. Besides, an 

estimation for can be found as the difference between Fn' (X) and ^^(A). For example, we can refer to the USE 
for the eigenvalues of the Hamiltonian calculated for various physical systems in "pure" states on the basis of the OM 

0-0... 

Considerable advantage of the USE for various applications in comparison to the asymptotic expansions is the 
possibility to investigate the qualitative peculiarities of the quantum system in the intermediate range of the parameters 
connected with Hamiltonian and the external conditions. At the same time, rather high precision of the USE in 
the zeroth order is of great use for practical calculations and it defines the rate of convergence of the consecutive 
approximations to the exact solution. 

Let us consider from this point of view the thermodynamic perturbation theory in the Schrodinger representation 
of the quantum statistics. Usually it is formulated for the free energy of the system |13) and the leading terms are 
the following 



n n m^n 

2 



1 



\ n / n 



V nn w n 



(3) 



Here we have introduced (3 = 1/kT where T is the temperature and k is the Boltzmann constant; F and Fo are the 

exact and approximate free energies respectively; E n °^ are the eigenvalues of the unperturbed Hamiltonian and V mn 
are the matrix elements of the perturbation operator with the following form of the total Hamiltonian 



H = Hq + fiV, 

and w n = exp [P(F - E n 0) )} is the unperturbed density matrix. 

Any partial sum of this series does not yield the free energy in the whole range of the temperature and the 
perturbation parameter /i even in the simplest cases. Let us illustrate this by means of the model Hamiltonian used 
earlier in for the convergence analysis of a usual perturbation series for "pure" states 
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H = H + nV = ~ (p 2 + x 2 ) + fix 2 ; 4°) =n+ 1 -- F = i In [2 sinh {(3/2)] ; 

V mn = (n + -)S m>n + - [y/(n + 1) (n + 2)<5 m) „ +2 + V"( ri ~ l)<^m,n-2], (4) 

where <5 mj „ is the Kroneker symbol. 

Certainly, the exact free energy is well-known for this model (F = A In [2 sinh ((3/2^/1 + 2//)]) but if we use these 
matrix elements in formula J3J), then rather simple calculation leads to the following result 



F^, (3) = i In [2 sinh 08/2)] + | coth | - 

^[1 + cothf + — !_(! + 2/3) + .•■]. (5) 
o ^ 2smh ^ 

It is evident that this series does not satisfy the USE criteria with respect to both essential parameters of the system. 
In the low temperature limit {(3 — > oo), formula JSJ) leads to the power series of fi for the ground state energy and this 
series diverges in the range of [i > 1/2 because of a singular point of the exact eigenvalue E n = yT+~2/7(n+ |) in the 
complex plane of fi 15] . On the other hand, Definition 1 breaks also for the free energy dependence on the parameter 
(3. When the temperature increases (f3 — ► 0), the second order correction has a stronger singularity (~ — /i 2 /4/3 2 ) than 
the exact free energy. 

Thus, our objective is to formulate another regular method, different from J2J|, which allows us finding the USE for 
the free energy of a quantum system with an arbitrary energy spectrum E n . 



III. THE OPERATOR METHOD FOR THE UNIFORMLY SUITABLE ESTIMATION OF THE 

EIGENVALUES 



As mentioned above, in order to obtain the USE for the partition function one has to solve two problems: 1) to find 
such an approximate representation for the energy levels of a system under consideration, which holds in the whole 
range of the Hamiltonian parameters, and 2) to perform approximate summation on quantum numbers so that the 
result remains uniformly suitable with any temperature. For the former purpose we shall use the operator method 
(OM) of approximate solution of Schrodinger equation. This method is described in details in review paper 0. Let 
us consider here only the basic expressions defining the consecutive approximations of the OM for the eigenvalues 
E n and eigenvectors \tp n ) of a Hamiltonian H of some quantum system (the index n may involve the whole set of 
quantum state numbers) 



H\ip n )=E n \ip n ). (6) 

Let us introduce the complete set of state vectors \n, u>), depending on the same quantum numbers and arbitrary 
parameters to. In the canonical perturbation theory (CPT) the eigenfunctions of some part Hq of the full Hamiltonian 
(zeroth order approximation Hamiltonian) are used as such a complete set. The Schrodinger equation for that part 
is assumed to be exactly solvable. In terms of the OM, the vectors |n, uj) are considered as a set of variational wave 
functions. Thus, the choice of these vectors is rather arbitrary and is defined by qualitative features of the system 
under consideration and by the possibility of rather simple calculation of the matrix elements of the total Hamiltonian 
H. Then we can represent the solution of equation © in the form of expansion 

\il> n ) = \n,u) + ^ f C nk \k,u;), (7) 
and choose the normalization of the exact solution as the following 



(^ n |n,o;) = l. (8) 

As a result, the initial eigenvalues problem © exactly reduces to the system of infinite number of nonlinear algebraic 
equations for E n and coefficients C n k 3] 
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E n — H nn + ^ C n kH n k] H n k{uj) — (n, u\ H \k, w) ; 

Cnm — \H 7nrn E 7l \ ^H mn -\- ^ y C nk H mk ^ , Tfl ^ 71. (9) 

As shown in [3|, the most effective method for solving this system of equations is an iterative scheme, in which the 
exact solutions are calculated as the limit of the sequences (and not by summation of series, as in the CPT) 

E n = lim £«; C nk = lim C%>, (10) 
and the consecutive approximations in (|10|) are calculated by means of recurrence relations 

k^n 

Cffl=0. (11) 

As shown for various physical systems Q - [U, recurrence relations (|llfl rapidly converge to the exact solutions, 
and the choice of the parameter uj influences only the rate of this convergence. In this paper we restrict ourselves 
with expressions for the energy after the second iteration (the first one coincides with the zeroth order approximation 
by definition) , as we are interested mainly not in strict numerical calculations, but in constructing of the USE for the 
eigenvalues, the partition function and the free energy. The second iteration leads to the following approximation 

E%\w) ~ £(°>(w) + AS„(w) - H nn - Wmm - H nn ]- l H nm H mn . (12) 

Strictly speaking, this expression is suitable if the diagonal matrix elements 

H kk {uj) ^ H nn (uj), 

for the considered value of the parameter to and for all states from the set \k,uj} in expansion l[7)l. 
If for some concrete values m and Lu m 

Hrnmi^ni) — E[ nn {iO rn ) , 

one should extract the state \m,tu m ) from expansion and solve the secular equation in the OM zeroth order 
approximation analogously to the canonical perturbation theory in case of degenerate states |l(j . Some applications 
of the OM for physical systems with such features were considered in papers 3 , @ • 

For systems with several degrees of freedom the state vector \ip n ,j) in equation may also depend on some 
additional index j. The energy level E n is supposed to be independent of this index. Actually it means that the wave 
function is also the eigenvector of some operator 



J\^n,j)=j\^n,j), [lH]=0, (13) 

which commutes with the Hamiltonian. As shown in paper in such a case the projection operator Tj for the 
state with fixed quantum number j should be used in series Q) in order to solve both equations © and ltT^|) 
simultaneously. It leads to unessential modification of recurrence equations l|ll|) but does not change the properties 
of the OM consecutive approximations. 

In this paper the generalization of the OM for quantum statistics is of main interest and thus we shall not consider 
the above mentioned more complicated systems. 
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In spite of the formal similarity of formula i|12|) and the expression obtained in the second order of the CPT (see 
formula @), the former one maintains some essential differences. First, the matrix elements of the operator H are 
calculated with the set of state vectors \n,uj) for the whole Hamiltonian (and not for its unperturbed part as in the 
CPT). According to this circumstance is essential for the convergence of the consequent approximations of the 
OM. In particular, this convergence depends on the numerical dimensionless parameter, defined by the ratio of the 
diagonal and nondiagonal matrix elements of the total Hamiltonian 



H,, 



k^n. (14) 



This value is independent of the physical parameters of the Hamiltonian and defines both the convergence of the 
consequent approximations of the OM, and the accuracy of the USE for the zeroth order approximation of the OM 
Q in accordance with Definition 1. 

Another specific feature of approximation (|12f> is that it depends on some undefined parameter u), which should be 
chosen so that to provide the best accuracy of the USE in the zeroth order approximation of the OM. It was shown 
in Q, that the best zeroth order approximation of the OM for the energy level with the quantum number n 

E n ~EW=EW(«> n ), (15) 
is achieved when uj — uj n is chosen as the solution of the equation 

^M=0. (16) 

OUJ 

It should be underlined, that equation l|l(jD is not the variational principle for the excited states. The point is that 
the artificial parameter uj defines actually the representation for the wave functions, and hence, the exact eigenvalues 
of the Hermitian Hamiltonian should not depend on its choice. The exact condition is 

^M=0. (17, 

Therefore, equation l|16fl can be considered as the OM zeroth order approximation for the exact condition (|17|l . 

It is essential that the optimal value of the parameter uo n depends on the quantum number of the considered state. 
This means that the orthonormal set of basic vectors used in expansion Q should be chosen for every state in different 
ways. As a result, the consecutive approximations of the OM (fill) have "local" character in the space of quantum 
numbers, i.e., they are calculated independently with different sets of state vectors for every state of the system. On 
the contrary, in terms of the CPT the corrections for all quantum numbers are defined by the same spectrum and 
basic state vectors of the unperturbed Hamiltonian. 

In this paper we use model physical systems in order to illustrate the applicability of the OM in statistical physics. 
In particular, let us consider the quantum anharmonic oscillator (QAO) with the Hamiltonian 

H = ^(p 2 + x 2 ) + f ix 2 + Xx 4 . (18) 

It is well known, that the CPT expansion for this system has zero radius of convergence with respect to the 
parameter A |l2^ . Therefore, the QAO is widely used for approbation of various nonperturbative methods, and all 
approximations can be compared with the detailed numerical analysis of this problem in [l^] ■ In the problem of the 
QAO the most adequate choice of the full set of state vectors for the OM algorithm is the set of eigenfunctions of the 
harmonic oscillator with an arbitrary frequency, playing a role of the parameter uj. At the same time, it is convenient 
to perform all the necessary calculations of the matrix elements in an algebraic form, using the representation of the 
secondary quantization 0] • In such a case, the variational parameter uj can be introduced directly into the Hamiltonian 
by means of creation a + and annihilation a operators 

(a + a + ), p = -ij~(a- a + ), 



2ZJ V V2 

a + a; [a,a + ]=l; n\n) — n\n) . (19) 
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Some fairly simple algebraic operations with expressions - (|16|) lead to the following results Q 



3A 



4 0) = ^>n + 1 + 2m)(2« + 1) + ^-(1 + 2n + 2n 2 )- 



6A(2n 2 + 2n + l) 

<4 - W „(l + 2/i) - 1 ^ + 1 ^ = 0; (20) 



AM 2 ) 



(n + l)(n + 2)[o; ra (l + 2/x - a; 2 ) + 2A(2n + 3)] s 
16w 2 [w„(w 2 + 1 + 2fi) + 3A(3 + 2n)] 
A 2 (n + l)(n+2)(n + 3)(n + 4) 



32w 2 K(w 2 + 1 + 2//) + 3A(5 + 2n)] ' 



AM 3 ) 



1 f 4 [(n + l)(n + 2)] 3/2 [2(2n + 3)A + u n {\ + 2/x) - t^] 3 
"256w 2 \ [3(2n + 3)A + cj„(l + 2^)+cj 3 ] 2 

2A [(n + l)(n + 2)] 3/2 ^(n + 3)(n + 4) [2(2n + 3)A + w n (l + 2/i) - w 3 



+ 



[3(2n + 3)A + co n (l + 2/x) + cj3] [3(2n + 5)A + w n (l + 2/*) + w 3 ] 
2A 2 (n + 3)(n + 4) [(n + l)(n + 2)] 3/2 [2(2n + 3)A + uj n {\ + 2/i) - ^ 3 ] 
[3(2n + 3)A + uj n (l + 2y) + lo^] [3(2n + 5)A + w n (l + 2/x) + w 3 ] 



| A 3 [(n + l)(n + 2)(n + 3)(n + 4)] 3/2 l 
[3(2n + 5)A + ^ n (l + 2^)+w 3 J 2 J 

It can be shown by strict comparison with numerical results |l2j | that simple analytical expressions (|20|l satisfy 
both definitions of the USE with £(°> ~ 0.03. Fig.l and Fig. 2 illustrate this statement. They show the results of 
calculation of the QAO energy levels in dependence on t he p arameter A and quantum number n, obtained by means 
of various methods: numerical calculations (scattered) [I^; the CPT; expansion in the strong coupling limit; the 
algebraic formulae (|2*U|> and (|2~T)> in the OM zeroth, second and third orders. The absolute values and relative errors 
of various calculation methods in dependence on the parameter A arc compared in Fig.l, and Fig. 2 shows the relative 
error in dependence on the quantum number n. One can see that only approximation for the functions E n {\) obtained 
in terms of the OM does satisfy the definition JTJ for the USE approximated energy levels of the QAO. 
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• Exact energy levels {Hioe et al, 1978) 
OM zeroth order 




T"| 1 1 1 | 1 1 1 | 1 1 1 1 — I I I I | 1 1 1 1 — I I I I | 



0,1 1 10 100 1000 

Anharmonicity parameter X 



Figure 1: Various estimations for the ground-state energy of the quantum anharmonic oscillator in dependence on the anhar- 
monicity parameter A. 
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Figure 2: Errors in calculation of the energy levels of the QAO in terms of various approximations considered as the functions 
of quantum numbers. 
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IV. OPERATOR REPRESENTATION FOR THE PARTITION FUNCTION AND THE CUMULANT 

EXPANSION 

Let us consider the partition function of some quantum system 



Z{(3) = ^g„exp[-/3£„ 



(22) 



n=0 



where E n is the energy level corresponding to a quantum number n and degeneracy g n . 

To obtain the USE for this function one has to solve two problems: 1) to construct the USE for the energy 
spectrum of a considered system that is valid for any quantum numbers, and 2) to obtain an approximate estimation 
for summation on quantum numbers suitable in the whole range of temperature. It should be pointed out, that for 
single-dimensional systems the latter problem is less important. It is sufficient to use expression l(T2"|) for the energy 
levels, obtained in terms of the OM, to derive the USE for thermodynamic values and to carry out direct numerical 
summation (see below) in (|22H . However, for systems with many degrees of freedom such summation is rather 
complicated numerical problem, therefore we shall discuss the possibility of constructing the USE for summation on 
quantum states. 

First, let us show that the partition function can be identically represented in the operator form as an average 
over some state vector. For that purpose one can use the basic set of the state vectors considered formally as the 
eigenfunctions of some excitation number operator 



h | n) = n | n) 



Then one can rewrite equation 1)22(1 as 



Z(f3) = (p*\exp[-pE(n) + 2(3* n - lnN(0*)] |/T 



(23) 



(24) 



Here |/3*) is the normalized "trial" state vector that depends on an arbitrary parameter /3* having the physical 
meaning of an effective inverse temperature for the equilibrium system of excitations considered in (|23|l . The actual 
value of (3* will be defined later from the condition of the best approximation for the partition function of a real 
system. By definition, vector |/3*) is of the following form 



= ^NJ2 v^exp[-/3*n] 

oo 

exp(-2/3*n) 



N 



.n=0 



(25) 



It is important that the state vector should not be considered as some expression for the mixed state cor- 
responding to Gibbs ensemble. Actually, equations ll'MI) - (|25H give the identical representation for the numerical 
(non-operator) value l|22[) in the form of quantum mechanical average using the auxiliary operator h. It is also nec- 
essary to stress that the operational representation for the partition function, analogous to (|2*4*| . can be obtained by 
choosing a more complicated trial vector instead of (|2*5|l . Such a vector should take into account specific features of 
the problem under consideration. In every particular case, such a choice may lead to a more accurate zeroth order 
approximation, but a simple vector (|25|l used in this paper allows us to construct the universal scheme for calculating 
the USE that is valid for an arbitrary quantum system. 

Besides, this representation is formulated without usage of some specific form for the degeneracy factors g n . In 
general case the exact calculation of this factor in advance could prove to be a rather complicated problem. But 
when the OM is used for the approximate calculation of the energy levels E n the factor g n can be calculated with the 
same accuracy by means of approximate solution of equation (|13J) which defines the degeneracy of the state with the 
quantum number n [3j. 

Now, we can use the cumulant expansion (CE) for calculating the average value l|24|l . Recall that the CE is valid 
for an arbitrary exponential operator when averaging over the normalized state vector |llj 



exp A 



exp 



(26) 
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where the cumulants K n are expressed in terms of the moments of the operator A. This expansion as a whole is an 
accurate one but at the same time every cumulant corresponds to the partial summation of a usual power series. As 
follows from , the first few terms in equation (j26|) are given by 



Kx = (A) ; K 



K a = (A i )-3(A)(A 



A 



2 A 



(27) 

The general procedure of the CE can be applied to the partition function in form l|24|l . Then successive approxi- 
mations for the partition function will be calculated taking into account the corresponding cumulants. In this paper, 
we restrict ourselves by the first and second cumulants which leads to the following approximation 



exp 



(A(n))-\nN+±((A 2 (n))-(A(n)) 2 ) 

A(n) = ~(3E{n) + 2/3*h 



(28) 



It is obvious that with any fixed number of cumulants in formula l|28|l the partition function Z(j3,/3*) depends 
on the artificial value 0* , which can be considered as a variational parameter, leading to the best approximation of 
the CE in some order. Particularly, the zeroth order approximation Zq([3,/3*) and the first correction Z\{J3, j3*) are 
defined by formulae 



Z = exp [-(3E(/3*) + 2/3* n - \nN((3*)] ; 



Z 1 = exp{^ P 2 {E 2 — E 2 ) — A(3(3*(En - En) + 4/3* a ( 



2( n 2 



}• 



(29) 



Here, all values are averaged over the "trial" distribution function corresponding to the ensemble of the excitations 
(12311 with the degeneracy of states g n and the effective inverse temperature f3* , as for example 



E(/3*) = N 9nE n exp[-20*n]; 

n=0 

N dN 



2 d(3* 



(30) 



Then the free energy calculated with the considered accuracy takes the following form 



F(J3) = -i In Z(J3) ~ F (0) + F^ + ■ ■ 
P 



= -~{\nZ [p*,(3] + In Zi\J3*,0\ + ■■■). 



(31) 



To obtain the estimation for the partition function of the given energy spectrum, the variational parameter j3* 
considered as the function of the real temperature (3 can be defined from the condition of the best approximation in 
the zeroth order of the CE. As follows from (|l(j|) . for this purpose one has to solve the equation 



8Z dE dn dN 



(32) 



One can use the energy spectrum, obtained in terms of the OM to calculate the moments E, E 2 with the distribution 
function l|3UI) . Then the totality of expressions (|12f> . (|16fl . (|29|l and (|32|l lead to the USE for the thermodynamic 
characteristics of a system with the accuracy to the terms of second order in approximation of the OM and the CE. 
In some cases, the additional approximation in these equations can also be applied. This approximation leads to the 
worst accuracy, but permits one to simplify significantly the evaluations. It is based on the following. 

Recall that condition l|16f> corresponds to the optimal choice of the parameter u n for the best approximation for the 
energy level with quantum number n in the zeroth order of the OM. But convergence of the successive approximations 



12 



of the OM exists for all values of u> |3j. At the same time, the partition function and the free energy are integral 
characteristics in the space of quantum numbers of a system. Thus, it is possible to expect that in order to construct 
the USE for these values it is sufficient to choose the parameter ui as the same value for all quantum numbers, and 
to use it for defining the best approximation condition for the free energy in the zeroth order of the OM and the CE. 
This means that for the calculation of the moments of the energy spectrum one has to take into account only the 
explicit dependence of the approximate energy spectrum on the quantum number n and disregard a more complicated 
dependence related to the choice of u> n from condition l|ltj|) • Hence, the values 



E(f3* , ^) ~ iV ]T g„ [E® H + AE n (w)] exp[-2/T n] ; 

oo 

W((3*,oj) ~ JV^ 5 „[i4 )H] 2 exp[-2/rn] (33) 

n=0 

depend on both the parameter (3* and the variational parameter u, which should be taken from the minimum condition 
for the free energy. As a result, instead of a system of algebraic equations l|29|) for u> n it is necessary to solve only one 
equation for the parameter lu 



— 9n [4 0) M] exp[-2/r n] = 0. (34) 

71=0 

The physical meaning of this equation is that the variational parameter u> should be chosen from the condition of 
the best approximation of the energy level with maximum occupancy corresponding to a given temperature. 

In this paper, we illustrate the efficiency of our method in obtaining the USE for some specific systems. For these 
systems, all the calculations can be performed analytically. In spite of the model pattern of the problems under 
consideration, they are rather widely used for investigations of thermodynamic properties of real molecular gases |14| . 
therefore, the set of analytical formulae considered below may be useful for nonperturbative analysis of anharmonic 
effects in such systems. 

It should be pointed out that in the case of harmonic oscillator (A = 0) in JTSJ, rather simple calculations show 
that, unlike the CPT (J3J), formulae l|33l) . 134|l lead to the exact expression for the partition function with any value of 
parameter ji. 

As mentioned above, the proposed method of calculating the USE for the partition function includes actually two 
components: the approximate summation over the quantum states and the approximate calculation of the energy 
spectrum. In order to illustrate how the first one works, let us obtain the USE for the partition function and the 
free energy for the system of quantum rotators. Such a system is used for simulation of the contribution of rotational 
degrees of freedom to the thermodynamic characteristics of diatomic molecular gases |13| . We can consider this system 
as a model with known energy spectrum and degeneracy, therefore it can be considered as a test for approbation of 
the cumulant expansion for approximate summation over the states with known spectrum. Hence, we should calculate 
the following partition function [13J 



oo 

Z r (f3,9 r ) = Z r (x) = J2( 2n + l)exp[-/36» r (n 2 + n)}. (35) 
n=0 

Here 8 r = h 2 /2I is the rotational temperature of the system measured in the energy units; / is the moment of 
inertia of the considered molecules. Let us introduce x as the dimensionless parameter, x = (39 r . 

One can find analytically the asymptotic formulae for small and large temperatures which are connected with the 
limit cases x 3> 1 and x <C 1 correspondingly 



Z r {x) ~ 1 + 3exp(-2a;) + . . . ; x > 1. (36) 
In the opposite limit case, wc use two terms of the Euler formula of summation 



oo 



f{y)dy- l -f{a)~^f{a), 
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which leads to 



Z r (x) 



f°° 11 
1 + Jo (2?/ + ^ ° Xp ^ x( - y2 + V ^ dy ~ 2 ~ 12 (2 ~ X) + ' ' ' = 



12' 



x < 1. 



(37) 



The presented expressions do not describe adequately the function Z r (x) in the range of intermediate values of the 
parameter x. Actually Z T (x) can be expressed through the Weierstrass function [24|, but it is rather complicated for 
analytical investigations. Let us show that the usage of the proposed method allows one to obtain the USE for this 
function in the entire range of x. 

In such a case all average values in formula H30J) can be directly expressed through moments of the "trial" distribution 
function 



N = 



£)(2n+l)e- 

.n=0 



-2/Tn 



(l-qf 



q = e- 2 ?* 



N^2n(2n+ l)e 

n=0 



i + q 

-2f3*n _ g(3 + g) , 



g(3 + 8g + g 2 ) 
(1 + qfl-qf 



(38) 



Using the presented formulae in l|29|) , we can obtain the following expression for the partition function in the zeroth 
order of the CE (the value < q < 1 can be considered as a variational parameter instead of (3* ) 



Zj?\x,q) = exp[<p(x,q)]; 



ip(x,q) 



6q 



3q 



{l-qf 



1-q 2 



■ In q — In 



(i-<z) 5 



This formula should be supplemented with an equation for calculating the function q 
condition (|32|) for the best choice of the zeroth order approximation 



(39) 

q(x), following from 



dtp 



(3 + 2 g + 3g 2 )(l-g) 
6(1 + q) 3 



\nq = 0. 



(40) 



3 + 2q + 3q 2 
(l-q)(l + qr 



-(3 + q)(l-q 2 ) 



In q — In 



1 + 9 



(41) 



System of equations 14UI) , 141fl should be considered as a parametric assignment of the function Z^ (x) , defining 
the zeroth order approximation for the partition function. 

The first order correction to the CE for the considered model spectrum is defined by the expression 



Z^(x,q) = exp[pi(a:,g)]; 



(n 2 + n) 2 — (n 2 + n) 



2x 



n(n 2 + n) — (n 2 + n)n 



\nq + 



In 2 q. 



Using equation 140J1 in the previous expression, one can find the following 



(42) 



Zi l \q) = exp[^ (<?)]; 

In 2 q. 



q 2 15 + 4q + 26q 2 + 4q 3 + 15g 4 , 2 



(l-q) 2 (l + qr 



(43) 



14 



The presented formulae define the USE for the partition function of the quantum rotator. In particular, they lead 
to the following results in the corresponding limit cases 

a;~-^, ip(q)~3q, (x) ~ 1 + 3 exp(-2a;) + . . . , x > 1; 

x~ (1 ~ g)2 , pfa) ~l + ln2-21n(l- ? ), Z<°)(x) ~ + 0.303 + ... , x « 1. (44) 

6 a; 

Comparing these results with corresponding asymptotic formulae H36(l . (|37|l . one can see, that the constructed 
approximation describes correctly the functional dependence of the exact partition function in both limit cases. At 
the same time in compliance with the common definition Q for the USE the accuracy of the estimation is defined by 
parameter £(°> ~ 0.1. The first-order correction (|43|l of the CE increases the accuracy of the estimation, in particular, 
when the value of x is small 



Z r {x) ~ Zl°\x)Zi 1 \x) ~ —— + 0.348+..., x<l. (45) 

x 

Fig. 3 represents the comparison of the results of evaluations in terms of exact and approximate analytical formulae 
in the zeroth and second orders of the CE. It demonstrates rather good accuracy of the zeroth order approximation of 
the CE and its uniform validity for all values of the dimensionless parameter, defining the properties of the discussed 
system. 




0,01 0,1 1 

Dimensionless Parameter x = p6 



Figure 3: Free energy of quantum rotator and error in its calculation in terms of the CE in zeroth and second order in 
dependence on the dimensionless parameter x. 
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V. THERMODYNAMICS OF THE QUANTUM ANHARMONIC OSCILLATOR 

The problem considered in this section illustrates the possibilities of the discussed method in the case when the 
approximate calculations are used for obtaining the USE for both the eigenvalues and summation over the states. 
For that purpose, we use the equilibrium statistical system of quantum anharmonic oscillators (QAO). This problem 
has been discussed in numerous papers concerning nonperturbative approximate methods (e.g. |l7j - |27|). but the 
approximation satisfying definitions and @ has not been obtained yet. As mentioned above, the solution of this 
problem can be useful for applications, in which one should take into account the anharmonicity of the molecular 
oscillations (for example, when considering either the thermodynamics of gases or the anharmonic effects in 
calculations of phonon spectra in crystals pjj). 

Let us consider the QAO with parameter fj, — 0. This can be done without loss of generality, since the parameter 
jj, can always be excluded from the Hamiltonian by the scaling transformation of the coordinate. We thus consider 



H=^(p 2 +x 2 )+\x 4 . 



(46) 



The USE for this system is defined by formulae 1(20(1 and 1(2111 . When obtaining the USE for the partition function 
and the free energy it is interesting to compare the results, which can be obtained by means of approximate calculations 
in two different ways: 1) direct numerical summation over the QAO states taking into account the OM energy levels 
and 2) evaluation of the USE in analytical form by means of the approximate summation over the states based on 
formulae 133(1 , 134(1 . To estimate the accuracy of our results we will use direct numerical calculation of the partition 
functi on 1221) (Z^ um ) for the QAO (g n = 1) taking into account 8 energy levels of the QAO, obtained numerically in 
paper [12j. It is essential to stress, that in real calculations with any fixed number of the QAO levels in the summation 
over the states, the numerical calculations lose their accuracy in the limit case of the high temperature (/3 — * 0), when 
the whole energy spectrum becomes essential (see left panel of Fig. 4). In these cases one can use for comparison the 
asymptotic expressions for the partition function which can be obtained analytically: 



Z A (0, A) = exp[-/3E n (A)] ~ 1 + e 



71=0 

1 - erf- 



3Ae~^ 
(1 -e~P) 2 



«i; 



{l-e-P) 2 



Z A {fi, A) ~ ]T exph/JA 1 / 3 ^], A » 1. 



(47) 



Here b n are known numerical coefficients of the asymptotic expansion for the QAO energy levels in the stronj 
coupling limit [12j . The analytical approximate expression for them in the OM zeroth order was obtained in paper 



3 4/3 l+2n + 2n 2 1/3 
""V [ (l + 2n) 2 1 • 



Fig. 4 represents the results of comparison of the free energy F 
expressions for it, which are defined as the following 



n u rn 
A 



(48) 

-1/ f3\n(Z^ um ) with various approximate 



^l 0) (/3,A) = -iln 



F| 1) (/3,A) = -iln 

id 01) (AA) = --Lm 



]Texp[-/3i4 >(A)] 

n=0 

£exph3(4°>+AE n )] 
n=0 

Y,eM-PE^](l-(3AE n ) 



.n=0 



(49) 



where the energy is calculated in different orders of the OM by formulae l|2U|) , 121(1 . We also compare our approxi- 
mations with the results of the canonical thermodynamic perturbation theory |13| and the direct summation over the 
exact eigenvalues for the first 8 levels of QAO calculated numerically in ^1 ■ 
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Figure 4: The free energy of the QAO using the OM zeroth and second orders and errors of these approximations in dependence 
on the anharmonicity parameter A (left panel) and on the inverse temperature j3 (right panel) in comparison with the results 
of the numerical calculations. Lines, corresponding to /3 = 1 on the upper part of the left panel, and all lines on the upper part 
of the right panel almost coincide. Results obtained in terms of the thermodynamic perturbation theory are also represented 
on the upper part of the left panel (they are not uniformly suitable and can not be shown on the right panel). 



The presented results show that the zeroth order approximation satisfies condition and defines the USE for 
the pointed values with a relative error ~ 0.1. If one takes into account the subsequent iteration of the OM, the 
accuracy of the estimation improves. At the same time, it is necessary to stress that the expansion of the exponent 
in a rather small value AE n used in expression (j3, A) breaks the USE conditions in the low temperature limit 
(/3 — > oo). This means that the improvement of the accuracy of the zeroth order approximation for the free energy in 
the subsequent iterations of the considered method has in general nonadditive character, unlike the thermodynamic 
CPT [13. 

Let us consider now the approximate expression for the partition function obtained in terms of the zeroth order 
approximation of both the OM and the CE due to formulae - ffiljl 

4° C) (/?,A)=exp[<K/?,/?*,A,w)]; 



<p(/3, (3* , A, w) = ~(3Ei 0) (/?* , w, A) + 2/3*n(/T ) - In N; 



2 



J V = ( 1 -,); B = rr A_;^=^iii p ; 9 = e-«- 
l,o .w_ 3A 



EL 0) = —(to 2 + l)(2n + 1) + — (1 + 2n + 2n 2 ). (50) 



17 



Variational parameters q and uj are defined from the minimum conditions for function ip, which lead to the following 
equations 



^ = o => -p " 2 + 1 + 3A(1 + g) +^ + 1 = 0- 

dq ' 2cj(1 — q) uj 2 (1 — q) 2 1 — q ' 

^=0,^^-^-6X^ = 0. (51) 

It is convenient not to solve numerically equations Q51JI. but to consider A(g, oS) and /3(q,w) as functions of variables 
q and w. Together with definition l|5U|l these functions assign in parametric form the functions Z^ C ^ (/?, A) and 

F^ C ^ ((3, A) = — 1//3 In i?^ ^ (/3, A) . The results of such a calculation and its comparison with the numerically calculated 
free energy are shown in Fig. 5. As one can see, the application of the CE to summation over the states keeps the 
zeroth order approximation of the method under consideration satisfying the criteria of the USE with a relative error 
~ 0.1. Fig.5 shows also that the values zf C) {f3,\) and F<£ C) (f3,\) are rather close to the corresponding values 

Z$\p,\) and F^\p,X) calculated by means of direct summation over the OM eigenvalues En (u) n ). In the latter 
case the parameter u n should be defined for each level separately. The former approach based on equations (33), (34) 
is much easier because the parameter uj is calculated only once for each temperature. It is especially important for 
systems with many degrees of freedom. 

To obtain the next approximation for the partition function Z A (/3, A) and the free energy F\ (/3, A) = 

-l/f3\nZ y A '(f3,\) it is necessary to take into account both the correction AE n for the energy levels obtained 
in terms of the OM, and the second cumulant in formula (|29|l for the approximate summation over the states. It is 
essential that in calculation of these corrections the same variational parameters q and ui are used, which were found 
in the zeroth order. Fig. 6 demonstrates the influence of these corrections on the accuracy of the USE. 
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Figure 5: Approximation for the free energy of the QAO using the CE zeroth order, direct numerical summation based on 
the OM zeroth order and exact results Jl2l together with the relative errors of such approximations in dependence on the 
anharmonicity parameter A for various values of inverse temperature j3 (left panel) and on the inverse temperature (3 for various 
values of the anharmonicity parameter A (right panel). Lines 4, 5 and 6 on the upper part of the left panel almost coincide. 
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Figure 6: Errors of various approximations for the free energy of QAO in dependence on the anharmonicity parameter A for 
various values of inverse temperature /3 (left panel) and on the inverse temperature j3 for various values of the anharmonicity 
parameter A (right panel). 

It is also of interest to consider the average values of observable physical characteristics of the system (e.g., energy) 
in dependence on the Hamiltonian parameters and temperature. For this purpose we use the numerical energy levels 
[l2|. the energy levels obtained in terms of the CPT and the OM in zeroth order. We consider also the expression 
for the average value of energy based on the cumulant expansion. Thus we obtain the following expressions for the 
average energy 



E(°) 



Z 



n 

exp[-/3i4 )] 



7 (0) 
J A n 



EiCPT) = ^_£4CPT) exp[ _^(CPT) ]; 



E( CE ) = - 



(CPT) 

_L d_ 

z(° c ) dp 



(52) 



The results of these calculations are presented in Fig. 7 which shows that our approach leads to the uniformly 
suitable estimation also for the average physical values unlike the canonical perturbation theory. As it was mentioned 
above, the noticeable deviation of the OM results from the numerical ones in the range of small (3 is explained by the 
fact that we could use only 8 eigenvalues found in ^1 f° r calculation of £X™ um ). It is also important to stress again 
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that the value calculated with the OM eigenvalues depending on the set of the parameters u„ is rather close to 
the value E^ CE ) depending on the only parameter uj from the equation (34). 




0,1 1 10 100 1000 0,1 1 10 

Anharmonicity Parameter X Inverse Temperature /? 



Figure 7: Average values of the energy of the QAO in dependence on the anharmonicity parameter A for various values of the 
inverse temperature ft (left panel) and on the inverse temperature [3 for various values of the anharmonicity parameter A (right 
panel) . 



VI. CONCLUSIONS 

In this paper, we have developed a nonperturbative method for calculation of the thermodynamic values of a 
quantum system. This has been achieved by combining the operator method of approximate solution of Schrodinger 
equation and the cumulant expansion for the summation over the quantum states. The method has been approved for 
the Boltzmann diatomic molecular gas in order to calculate the partition function and the free energy defined by the 
molecular internal degrees of freedom. We have used some realistic models for the molecular movement (a quantum 
rotator and a quantum anharmonic oscillator) and found the uniformly suitable estimation for the thermodynamic 
values. This estimation tends asymptotically to the exact expansions in limit cases of both temperature and the 
Hamiltonian parameters. Besides, the zeroth order approximation of the proposed method is in close agreement 
with the exact results (with the relative error no more than 0.1) in the whole range of both temperature and the 
Hamiltonian parameters. A systematic procedure for calculation of the subsequent corrections has been formulated 
and the second order corrections have been found to improve the accuracy of the estimation. It has been also 
shown that application of the cumulant expansion for the summation over the system states permits one to calculate 
directly the thermodynamic values without preliminary high precision estimation of the whole set of the Hamiltonian 
eigenvalues. It is especially important for application of the proposed algorithm for the systems with many degrees 



21 



of freedom. 

Certainly, the strict proof of convergence of the formulated method is of special interest. Unfortunately it seems to 
be a very complicated mathematical problem in general case. But approbation of the method for a series of model 
systems in this paper can be considered as a qualitative argument for possibility of applying of this approach to real 
physical problems which we suppose to analyse in the forthcoming papers. 
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VII. CAPTIONS FOR FIGURES 

Figure 1. Various estimations for the ground-state energy of the quantum anharmonic oscillator in dependence on 
the anharmonicity parameter A. 

Figure 2. Errors in calculation of the energy levels of the QAO in terms of various approximations considered as 
the functions of quantum numbers. 

Figure 3. Free energy of quantum rotator and error in its calculation in terms of the CE in zeroth and second order 
in dependence on the dimensionless parameter x. 

Figure 4. The free energy of the QAO using the OM zeroth and second orders and errors of these approximations 
in dependence on the anharmonicity parameter A (left panel) and on the inverse temperature (3 (right panel) in 
comparison with the results of the numerical calculations. Lines, corresponding to (3 = 1 on the upper part of the 
left panel, and all lines on the upper part of the right panel almost coincide. Results obtained in terms of the 
thermodynamic perturbation theory are also represented on the upper part of the left panel (they are not uniformly 
suitable and can not be shown on the right panel). 

Figure 5. Approximation for the free energy of the QAO using the CE zeroth order, direct numerical summation 
based on the OM zeroth order and exact results 01 together with the relative errors of such approximations in 
dependence on the anharmonicity parameter A for various values of inverse temperature (3 (left panel) and on the 
inverse temperature (3 for various values of the anharmonicity parameter A (right panel). Lines 4, 5 and 6 on the 
upper part of the left panel almost coincide. 

Figure 6. Errors of various approximations for the free energy of QAO in dependence on the anharmonicity 
parameter A for various values of inverse temperature /3 (left panel) and on the inverse temperature (3 for various 
values of the anharmonicity parameter A (right panel). 

Figure 7. Average values of the energy of the QAO in dependence on the anharmonicity parameter A for various 
values of inverse temperature (3 (left panel) and on the inverse temperature (3 for various values of the anharmonicity 
parameter A (right panel). 



